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Abstract. A new path-integral theory is developed to calculate the photoemission spectra 
(PES) of correlated many-electron systems. The application to the study on Bi2Sr2CaCu208 
(Bi2212) and boron-doped diamond (BDD) is discussed in details. It is found that the 
isotopic shift in the angle-resolved photoemission spectra of Bi2212 is due to the off-diagonal 
quadratic electron-phonon (e-ph) coupling, whereas the presence of electron-electron repulsion 
partially suppresses this effect. For the BDD, a semiconductor-metal phase transition, which is 
induced by increasing the e-ph coupling and dopant concentration, is reproduced by our theory. 
Additionally, the presence of Fermi edge and phonon step-like structure in PES is found to be 
due to a co-existence of itinerant and localized electronic states in BDD. 



1. Introduction 

The role of electron-phonon (e-ph) interaction in the high-T c superconductivity has received 
considerable attention since the discovery of oxygen isotope effect of E^S^CaC^Og (Bi2212) 
in the angle-resolved photoemission spectra (ARPES)[1] [2]. The experimental data show that 
electronic band is shifted slightly with the 16 0/ 18 substitution, signifying the existence of e-ph 
interplay in cuprates. Besides, theoretically clarifying this effect is of great significance, for the 
energy scale of shifts reported by the two groups in Refs. pQ and [2] seem quite inconsistent 
with each other, and many questions still remain up to now. In order to have an insight 
into the isotope effect, in this work, we develop a new path- integral theory to calculate the 
photoemission spectra (PES) of cuprate superconductors, in which the electron-electron (e-e) 
and e-ph correlations are treated on an equal footing. This theory is also applicable to other 
kind correlated materials. As an example, here, we also study the PES of boron-doped diamond 
(BDD), which undertakes a semiconductor-metal phase transition on doping, and becomes a 
superconductor with the temperature decreases [3]. The details of our theory will be presented 
in the next section, in connection with the study on isotope effect. Calculation and discussion 
on PES of BDD are included in Section 3. A concluding remark can be found in the Summary. 



2. Isotopic shift in ARPES of Bi 2 Sr 2 CaCu20 8 

2.1. Model for CuOi plane of cuprate superconductor 

In the Cu02 plane of cuprates, the electronic transfer is modulated by the vibration of oxygen 
atoms between the initial and final Cu sites (see in Fig. 1), resulting in an off-diagonal type 



e-ph coupling. In order to qualitatively clarify the isotope effect of Bi2212, we start from a 
half-filled Hamiltonian including the e-e repulsion and the above mentioned off-diagonal e-ph 
coupling (ft, = 1 and ks = 1 throughout this paper): 



where aj (a/ CT ) is the creation (annihilation) operator of an electron with spin a at the Cu site 
I on a square lattice (Fig. 1). The electrons hop between two nearest neighboring Cu sites, 
denoted by (I, I'), with a transfer energy t(l, I'). U is the strength of Coulomb repulsion between 
two electrons on the same Cu site with opposite spins. The oxygen phonon is assumed to be of 
the Einstein type with a frequency loq and a mass m. A (= 1 + Am/m) is the mass change factor 
of phonon due to the isotope substitution. In the third term, qui is the dimensionless coordinate 
operator of the oxygen phonon locating between the nearest-neighboring Cu sites I and I', and 
the sum denoted by {I, I') just means a summation over all the phonon sites in the lattice. 



Figure 1. Lattice structure of Cu02 conduction 
plane in cuprates. The copper atom (white 
circle) is on the l-th site of a simple square 
lattice, where the electrons (arrows) reside. The 
oxygen atom (black circle) is located between the 
nearest-neighboring Cu sites, and qu> denotes its 
displacement from the equilibrium position. 



^ 3 q ^ 4 

o • o 



1 q 12 

" 1 r 



o • o 

Cu O 



In the conduction plane of Cu02, the electronic hopping integral t(l,l') can be expanded to 
the second order terms with respect to the phonon displacements qw as 

t(l,l')=t + sqf l ,, (2) 

where to is the bare hopping energy and s the off-diagonal quadratic e-ph coupling constant. 
Here we note the linear e-ph coupling does not occur owing to the lattice symmetry of present 
model. Whereas the inter-site e-e interaction is included in the screened values of to an d s. 

2.2. Path-integral Monte Carlo method 

In this section, we develop a path-integral theory for a model with both e-e and e-ph interactions. 
By making use of the Trotter's decoupling formula, the Boltzmann operator is written as, 



lim (e- AH ---e' AH ). (3) 



Applying the Hubbard-Stratonovitch transformation [3] and the Gaussian integral formula[5j, 
we can decouple the two-body parts, so that the e-e and e-ph correlated terms are replaced by 
a two-fold summation over the auxiliary spin and lattice configurations, which is the so-called 
path-integral. In this way, the Boltzmann operator is rewritten into the path-integral form as, 

e~ PH -» y"Px(r + exp|-jfV[/i(T,z^ (4) 
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(5) 
(6) 



Here, x m and x q correspond to the auxiliary spin and lattice field, respectively, J T>x symbolically 
denotes the integrals over the path x synthesized by x m and x q , and \x q ) is the eigenstate of 
phonon. A is the time interval of the Trotter's formula, (3 = 1/T, and T is the absolute 
temperature. T + in Eq. (4) is the time ordering operator. 

Then the time evolution operator [= R(t, x)\ along a path x is defined as 



R(t, x) = T + exp 



dT'h(T',X m ,X q ) 



(7) 



In terms of the Boltzmann operator (4) and time evolution operator (7), we define the free 
energy [= 3>(cc)] of the given path as 



While, the partition function (= Z) and total free energy (= are given as 



(8) 
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According to Refs. [1] and [5], we also define the one-body Green's function [= G a (lr, I't', x)] 
on a path x as 

G a (lT,l'T',x) = -sign(r - T')(T + a la (T)a\ l(T (T')) x , (10) 
where ai a {r) is the Heisenberg representation of a\ a . It is really time-dependent and defined by 

ai a (r) = R~ x (t, x)ai a R(r, x). (11) 
Meanwhile, the ordinary Green's function [= G a (l,r)] can be obtained by the path-integral as 

G a {l -/',r-r / ) = || Vxe-f>*V>G a {lT, I't', x). (12) 

This path-integral is evaluated by the quantum Monte Carlo (QMC) simulation method. 

If the QMC data of Green's function G a (l,r) is obtained, we can immediately calculate its 
Fourier component [= G a (k,r)] as 



1 



ik-R 



(13) 



where k is the momentum of the outgoing photo-electron. From this Fourier component G a (k, r) , 
we derive the momentum-specified spectral function [= j4 CT (k, u)] by solving the integral equation 
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(14) 



2.3. Isotope substitution induced band shift in ARPES 

We now present the QMC results on a 4x4 square lattice, where to is set as the unit of 
energy, and o>o=1.0 is used. For the QMC simulation, we impose a little large isotopic mass 
enhancement, Ao=l and A=2, to suppress the numerical error. In this calculation, we determine 
the binding energy ek by the moment analysis of the spectral function as ek = I dwA a {\a, lu)lu. 
Correspondingly, the isotope induced band shift is calculated by Aek = ek(A) — ek(Ao). 




Figure 2. (a) The variation of Aer/Ao; with s on a 4x4 square lattice, when [7=2.0, 
(3=10, Ao=l, A=2. The filled circles are from QMC, and the solid curve from HFA as a 
guide for eyes, (b) The variations of Aer/Ao; with [7 on a 4x4 square lattice at 0=10, 
Ao=l and A=2. Three different values of s are used to show the s-dependence of Aer/Ao;. 
The discrete symbols are the results of QMC, and continuous curves by HFA as a reference. 



In Fig. 2(a), we plot the ratio Aer/Ao; versus s, at [7=2.0 and (3=10, where Aer is the 
band shift at the T point of Brillouin zone [kr=(0,0)], and Ao; is the isotopic change of phonon 
energy. The filled circles are calculated by QMC, and the solid curve by the mean-field theory 
with Hartree-Fork approximation (HFA) as a guide for eyes. Here both theories figure out an 
increase of Aer/Ao; with s, which means if the e-ph coupling is strong enough, a large band 
shift can be generated in a small cost of Ao;. In Fig. 2(b), the ratio Aer/Ao; versus U are shown 
for three different s's, where the discrete symbols and continuous curves are the QMC and HFA 
results, respectively. One can see the ratio Aer/Ao; increases with s. Meanwhile, for a fixed 
s, the ratio declines slightly as U increases, indicating that the band shift is owing to the e-ph 
coupling, whereas the presence of U partially reduces this effect. According to Fig. 2, the band 
shift thus can be regarded as a measure of the effective e-ph coupling strength in the system. 

3. Co-existence of localization and itineracy of electrons in boron-doped diamond 

It is well known that the pristine diamond is a big band gap insulator. Lightly doped with boron, 
it shows ap-type character with an activation energy about 0.37 eV[6J. Recently, the research on 
BDD has become highly attractive since the remarkable discovery of superconductivity in this 
material[3]. Accompanied with the superconducting phase, a semiconductor-metal transition 
also occurs when the doping percentage is increased to certain level. In the normal metallic 
state, K. Ishizaka et al.[7\ declared the observation of a step-like multi-phonon satellite structure 
in the valence band PES, approximately distributed periodically in 0.150 eV below the Fermi 
level, in addition to the emergence of a clear Fermi edge. This periodic structure in PES reminds 
us of a similarity to the case of localized electron [8], wherein the coupling between electron and 
Einstein phonons characterizes the spectra with discrete peaks at equal distance. Moreover, 
the Fermi edge and the step-like structure are observed together, probably originating from a 
co-existence of the two basic properties of electrons: itineracy and localization. In order to 
clarify this so-called co-existence theoretically, we apply our path-integral theory to the many- 
impurity-Holstein (MIH) model on a doped simple cubic lattice to derive its spectral density. 



3.1. Model and methods 

The MIH model includes the following two properties. One is the disorder of the system, 
that some atoms are replaced by dopant ones in a certain ratio. The other is the coupling 
between electrons and Einstein phonons, being the simplest description of e-ph interactions. Its 
Hamiltonian is given as, 




-S^2Qi(ni a -ni/2), n la = a\ a a ia , a = a or /3, (15) 

where t is the transfer energy. a[ and a\ a are the creation and annihilation operators of electron 
with spin a at site I. Electrons can hop only between the nearest neighboring sites expressed by 
(I, I'). fM stands for the chemical potential of electrons and A e is the potential difference between 
after and before substitution at the doped sites labeled by Iq. Qi is the dimensionless coordinate 
operator for the phonon at site I with frequency ujq. S denotes the e-ph coupling constant, ni 
is the average electron number at site I. To simplify the problem, we just consider the coupling 
at the doped sites hereafter, because in the pure diamond there is no evidence that the satellite 
structure appears, which suggests the coupling is important only after doping. 

By applying our path-integral theory to the disordered system, and following the formulation 
in Section 2.2, we obtain an electronic Green's function similar to Eq. 12, 

G a (l, r) = ^J Vxe-^G a (lT, 10, x). (16) 

In the numerical calculation, the path-integral of Green's function in Eq. 16 is also performed by 
the QMC simulation method as before. After averaging this site-dependent Green's function over 
all the N sites of the system, the photoemission spectral function [= N a (u)] can be reproduced 
through the analytic continuation as, 

-XG.M—f^j^N.iuVu. (17) 

Finally, after imposing the Fermi-Dirac function f{oj) = l/[exp(/?u) + 1], the PES intensity is 
given as = ^2 a N a (uj)f(uj), which can be compared with the experimental data. 

3.2. Results and discussions 

Since we focus on the spectral region close to the Fermi level, we use a simple cubic lattice of 
4 x 4 x 4 in real calculation without paying much attention to the detail of carbon valence band. 
As we just count the e-ph coupling at the doped sites, the phonon effect is not obvious in the 
whole system PES after averaging over all sites. For this reason, in the following, we will also 
present the PES of doped sites to illustrate the phonon effect. 

Fig. 3 shows the spectra for the simple cubic lattice at light (3.12%) and heavy (12.5%) 
doping rates, with a weak e-ph coupling of »S=0.25 eV. From the whole system PES (full lines), 
one can clearly see the emergence of a Fermi edge with the increase of doping percentage. In 
the lightly doped case, the impurity levels are a little above the top of valence band, indicating 
the system is a semiconductor. While in the heavily doped case, the impurity band expands to 
overlap the top of valence band and closes the activation gap, leading to a semiconductor-metal 
phase transition. In the spectra of boron-doped sites (dotted lines), one notices that, in addition 
to a broad hump due to the so-called Fano effect [9], a satellite structure also appears in each case 
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Figure 3. The QMC-calculated spectra for 
a simple cubic lattice of light (left) and heavy 
(right) doping rates, at a weak e-ph coupling. 
Zero is the position of Fermi level. Full lines 
denote the PES of whole system, and dotted 
lines the spectra of doped sites. 
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Figure 4. The co-existence of Fermi edge 
and step- like structure in the PES, at a low 
doping rate and a strong e-ph coupling, by 
QMC. Full (dotted) line is the spectrum of 
whole system (doped sites). The inset zooms 
in the solid line near the Fermi level. 



a little below the Fermi level. We should note this side structure is just the quantum phonon 
peak due to the e-ph coupling. Meanwhile, because of the Fano effect, this phonon structure is 
modified by the impurity band to have a shoulder-like shape, rather than a (5-like form as in the 
localized electron case of Ref. [8]. 

In Fig. 4, we show the spectra of a system with a strong e-ph coupling of 5=0.50 eV, and 
a low doping rate of 3.12%. Comparing with Fig. 3(a), one clearly finds that, a large e-ph 
coupling at the doped sites greatly contributes to the expansion of impurity band, so that the 
Fermi edge is clearly observed in the whole system PES (full line). Moreover, the second phonon 
side peak also appears due to the multi-phonon processes aroused by the strong e-ph coupling, 
corresponding to a step-like structure in the whole system PES, as affirmed in the inset. 

4. Summary 

We develop a new path-integral theory to calculate the PES of correlated many-electron systems. 
The isotopic shift in the ARPES of Bi2212 is investigated by this theory based on a model 
including both e-e and off-diagonal quadratic e-ph interactions. Our calculation demonstrates 
that the band shift is primarily triggered by the e-ph coupling, while the presence of e-e 
repulsion tends to suppress this effect. We also apply this theory to the MIH model on a 
doped cubic lattice to clarify the spectral properties of BDD. It is clearly shown in the PES that 
a semiconductor-metal phase transition takes place due to the increases of e-ph coupling and 
dopant concentration. Furthermore, the presence of Fermi edge and phonon step-like structure 
indicates the co-existence of two basic characters of electron, itineracy and localization, in BDD. 
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